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The  primary  objectives  of  this  research  were  to  develop 
capabilities  for  estimating  the  degree  of  revegetation 
success  to  be  expected  under  a  wide  variety  of  climatic 
conditions,  soil  and  spoil  properties,  and  revegetation 
treatments  utilizing  site-specific  revegetation  data  and 
information  from  most  of  the  coal  surface  mines  in  the 
interior  West.  The  more  important  study  developments 
are  these:  A  strong  conceptual  framework  for  evaluating 
the  success  of  proposed  vegetative  rehabilitation  efforts 
on  areas  to  be  surface  mined.  Site-specific  maps  to  pro- 
vide precipitation,  growing  season  length,  soil,  and  vege- 
tation type  information,  critical  to  the  evaluation  system. 
Results  of  most  existing  surface-mine  rehabilitation 
efforts  in  the  interior  West  through  1976  were  surveyed 
as  a  basis  for  evaluating  revegetation  success.  Interim 
predictors  (models)  were  developed,  based  on  these 
survey  data,  for  forage  production  and  cover  density 
potentials  on  undisturbed  sites  adjacent  to  the  study 
mines  and  mined  areas  with  various  combinations  of 
planting  and  postplanting  methods  and  treatments 
designed  to  enhance  the  total  rehabilitation  effort  for 
native  species  and  introduced  species. 

This  research  shows  that  practical  criteria  for  measur- 
ing revegetation  success,  namely,  the  amount  of  forage 
produced  and  the  density  of  plant  cover  developed,  are 
affected  significantly  by  at  least  two  major  climatic 
factors  (precipitation  and  growing  season)  that  are  not 
readily  susceptible  to  alteration;  by  three  properties  of 
spoil  materials  (potassium,  sodium,  and  pH)  that  are 
subject  to  limited  modification  through  management. 
Revegetation  success  is  also  influenced  by  seven  revege- 
tation treatments,  each  of  which  provides  at  least  two 
management  alternatives,  and  by  the  age  of  the  vegeta- 
tion. These  characteristics  account  for  about  one-half  to 
three-fourths  of  the  total  variance  in  forage  production 
and  plant  cover  density  in  the  prediction  models. 
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INTRODUCTION 

An  emerging  problem  in  the  interior  West  is  the  ad- 
verse effect  on  environmental  quality  of  spoils  left  in  the 
wake  of  surface  mining  for  coal  (cover  photo,  "Mining"). 
Needed  are  criteria  and  guides  for  predicting  revegetation 
potentials  on  various  kinds  of  surface-mined  land.  Equally 
important  is  the  need  to  define  and  prescribe  revegetation 
treatments,  as  well  as  posttreatment  management  mea- 
sures, for  such  land. 

Passage  and  implementation  of  the  Surface  Mining 
Control  and  Reclamation  Act  (Public  Law  95-87)  and 
attendant  regulations  have  placed  a  new  emphasis  on 
revegetation  of  spoil  materials  from  coal  surface  mines 
in  North  Dakota,  Montana,  Wyoming,  Colorado,  Utah, 
New  Mexico,  and  Arizona.  Although  most  coal  com- 
panies in  the  interior  West  are  applying  revegetation 
practices  to  surface  mines,  requirements  of  the  new  Fed- 
eral Reclamation  Act,  in  many  cases,  will  require  a 
reassessment  and  evaluation  of  techniques  and 
methodologies  employed  in  present  revegetation  activi- 
ties. Currently,  many  research  activities  are  underway  to 
determine  the  best  "mix"  of  cultural  practices  and  plant 
species  needed  to  satisfactorily  revegetate  disturbed 
land  (cover  photo,"Revegetation").  In  the  interior  West  not 
enough  time  has  elapsed,  however,  since  surface  mine  re- 
vegetation activities  began  to  give  assurance  that  any  par- 
ticular combination  of  revegetation  methods  will  be  suc- 
cessful in  the  long  run.  Consequently,  mining  applicants, 
as  well  as  administrators  granting  approval  to  mine,  can 
only  guess  on  the  basis  of  limited  experience  whether  the 
required  reclamation  standards  can  be  met  or  not. 

In  view  of  this  uncertainty  regarding  the  probability 
for  successfully  revegetating  western  surface  coal  mines, 
an  investigation,  financed  jointly  by  the  Forest  Service's 
Surface  Environment  and  Mining  (SEAM)  project,  the 


Environmental  Protection  Agency,  and  the  Fish  and  Wild- 
life Service,  was  begun  in  1976  to  identify  criteria  for 
measuring  the  success  of  revegetation,  to  evaluate  past 
and  ongoing  revegetation  efforts  at  most  of  the  major 
surface  coal  mines  in  the  interior  West,  and  to  develop  a 
capability  for  predicting  the  probable  degree  of  revegeta- 
tion success  expected  on  coal  lands  that  are  surface 
mined  in  the  future.  That  investigation  and  its  results  are 
the  subject  of  this  paper. 

INVESTIGATIVE  METHODS 

Criteria  as  set  forth  by  the  Surface  Mining  Control  and 
Reclamation  Act  (P.L  95-87,  Section  515  (6)(19))  require 
coal  mine  operators  to  "establish  ...  a  diverse,  effective, 
and  permanent  vegetative  cover  of  the  same  seasonal 
variety  native  to  the  area  of  land  to  be  affected  and 
capable  of  self-regeneration  and  plant  succession  at 
least  equal  in  extent  of  cover  to  the  native  vegetation  of 
the  area;  except  that  introduced  species  may  be  used  in 
the  revegetation  process  where  desirable  and  necessary 
to  achieve  the  approved  post-mining  land  use  plan." 

Accordingly,  investigative  methods  were  based  on  the 
assumption  that,  unless  some  other  revegetation  objec- 
tive is  defined,  the  primary  goal  of  revegetation  on 
surface-mined  coal  areas  in  the  interior  West  is  to  estab- 
lish a  productive  and  protective  cover  of  durable  plants, 
consisting  primarily  of  species  adapted  to  and  character- 
istic of  similar,  but  unmined  areas.  It  was  further 
assumed  that  differences  in  the  degree  of  revegetation 
success  to  date  on  surface-mined  coal  lands  in  the  West 
should  be  related  to  variations  in  natural  climatic  com- 
ponents, changes  in  site-specific  physical  and  biological 
characteristics,  and  differences  in  the  revegetation 
methods  used. 
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The  degree  to  which  plant  cover  is  established,  either 
in  natural  or  revegetated  stands,  can  be  measured  and 
evaluated  in  a  nunnber  of  ways.  One  of  the  most  impor- 
tant measures  of  success  in  plant  cover  establishment  is 
the  capability  of  the  vegetation  to  produce  aboveground 
biomass  for  forage  or  some  other  useful  purpose. 
Another  important  measure  of  success  is  the  capability 
of  the  vegetation  to  produce  ground  cover  for  protection 
of  the  soil  against  the  erosive  forces  of  raindrops  and 
surface  runoff.  Accordingly,  the  degree  of  success  of 
vegetation  reestablishment  was  measured  in  terms  of 
total  weight  of  aboveground  biomass  and  total  percent 
density  of  ground  cover  consisting  of  plant  basal  area 
and  accumulated  litter.  Both  of  these  measures  are  com- 
patible with  the  legislative  constraints  under  which  mine 
operators  must  evaluate  reclamation  success. 

With  these  constraints  in  mind  and  exceedingly  strong 
agency  pressures  for  at  least  an  interim,  but  "immedi- 
ate," process  for  evaluating  proposed  revegetation 
actions  on  areas  to  be  surface  mined,  a  study  of  existing 
mine  revegetation  efforts  in  the  interior  West  was  under- 
taken. The  planned  strategy  was  to  provide  information 
on  forage  production  and  cover  for  unmined  areas  to 
define  undisturbed  site  potentials,  and  to  provide  similar 
information  for  comparable  surface-mined  sites  to 
measure  the  effectiveness  of  vegetative  rehabilitation 
treatments  applied  to  date.  Those  combinations  of  treat- 
ments applied  to  mined  sites  and  producing  forage  and 
cover  conditions  at  least  equivalent  to  the  potentials  on 
associated  undisturbed  sites  would  then  be  considered 
satisfactory  for  the  legislative  objectives  above. 

During  the  growing  seasons  of  1976  and  1977,  data 
were  obtained  from  28  of  the  34  major  coal  surface  mines 
then  located  throughout  the  surface-minable  coal  areas 
of  the  West.  These  data  provide  information  that  de- 
scribes important  climatic  features,  physical  and  biologi- 
cal characteristics  of  each  site,  treatment  measures 
employed  to  effect  revegetation,  the  age  of  the  planting 
on  each  revegetated  area,  and  amounts  of  forage  and 
ground  cover  density  developed  by  both  native  and  intro- 
duced types  of  vegetation  on  these  areas.  Similar  infor- 
mation, except  for  the  age  of  the  vegetation,  was  also  ob- 
tained for  unmined  areas  near  each  mine  that  had  been 
long-undisturbed  and  were  characterized  by  predomi- 
nantly native  vegetation. 

Information  about  general  climatic  features  was 
obtained  from  tabulations  of  State  climatic  data  and  con- 
sisted of  total  annual  precipitation,  growing  season  pre- 
cipitation, and  length  of  the  frost-free  growing  season. 
Climatic  information  of  significance  appears  on  the 
maps  of  appendix  B. 

Data  concerned  with  site-specific  physical  and  bio- 
logical factors  were  obtained  from  176  100-foot  transects 
selected  to  be  representative  of  the  mined  and  adjacent 
unmined  study  areas.  Measurements  included  the 
aspect,  slope  steepness,  and  elevation  of  each  transect. 
Soil  samples  from  the  unmined  areas  and  spoil  samples 
from  the  mined  and  revegetated  areas  were  obtained  to  a 
depth  of  8  inches  (20  cm)  below  the  ground  surface  and 
were  analyzed  for  texture,  conductivity,  nitrate  nitrogen, 
phosphorus,  potassium,  sodium,  calcium,  magnesium, 
acidity  (pH),  sodium  adsorption  ratio,  and  saturation  per- 
centage. 


Along  each  study  transect,  fifty  1-ft^  (0.1  m^)  plots  were 
located  randomly  for  use  in  determining  species  com- 
position, aboveground  biomass  or  forage  production,  and 
plant  cover  density.  All  current  growth  of  perennial  plants 
within  these  plots  was  clipped  to  a  height  of  one-half 
inch  (1.27  cm)  above  ground,  bagged  by  species,  oven- 
dried  for  24  hours,  and  weighed. 

The  age  of  each  revegetated  area,  expressed  as 
number  of  years  since  it  was  planted,  and  information 
concerning  the  treatments  applied  during  and  subse- 
quent to  its  establishment  were  obtained  through  consul- 
tation with  reclamation  personnel  employed  at  each  of 
the  mines.  These  treatments  reflected  differences  in 
tillage,  seeding  methods,  topsoiling,  fertilizing,  supple- 
mental irrigation,  mulching,  and  time  of  seeding. 

All  of  the  data  obtained  from  176  transects  on  the  28 
mines  studied  made  up  the  information  base  available 
for  analyses. 

DATA  AND  INFORMATION  ANALYSES 

Six  series  of  analyses  were  made.  These  analyses 
followed  multiple  regression  strategies  for  estimating 
forage  production  and  plant  cover  density  as  functions  of 
climatic  components,  growing  medium  characteristics, 
and  revegetation  treatment  alternatives.  The  simple 
linear  effects  of  independent  variables  on  forage  pro- 
duction and  cover  density  were  screened  statistically  in 
all  combinations  as  a  means  of  isolating  the  stronger 
variables  for  use  in  synthesizing  final  models. 

Forage  Production  Model  for  Unmined  Areas 

Eighty-three  of  the  176  transects  sampled  during  this 
investigation  were  located  on  unmined  and  otherwise 
long-undisturbed  areas  adjacent  to  each  of  the  28  mines 
studied.  For  these  areas,  annual  precipitation  amount 
and  growing  season  length  each  added  significantly 
(Pr  <  0.005)  to  the  regression  for  forage  production.  Soil 
potassium  content,  while  not  adding  significantly  to  the 
regression  after  fitting  precipitation  and  growing  sea- 
sons, did  display  unusual  strength  in  the  short-to-medium 
length  growing  season  and  medium-to-high  precipitation 
range;  so  it  was  retained  in  the  model. 

Strong  interactions,  not  likely  to  be  well  represented  by 
the  simple  linear  additive  effects  initially  screened,  were 
expected  to  exist  among  these  variables.  Accordingly, 
attention  was  focused  on  the  interactive  effects  of  these 
three  variables  on  forage  production.  Forage  production 
was  expected  to  increase  upward  concavely  with  increas- 
ing precipitation,  to  reach  a  peak  somewhere  within  the 
broad  range  of  growing  season  lengths  encountered,  and 
to  increase  with  increasing  amounts  of  potassium  in  the 
soil.  The  interactive  effects  of  these  variables  were 
modeled  under  the  constraints  of  expectation,  following 
Jensen  and  Homeyer  (1970,  1971)  and  Jensen  (1973,  1976, 
1979).  More  expressly,  expected  trends  were  adjusted  to 
the  partitioned  data  graphically,  coordinated  to  arrive  at 
the  interactive  relation,  and,  simultaneously,  formulated 
mathematically.  The  resulting  model  was  adjusted  to  the 
unmined-transect  source  data  (n  =  83)  by  least  squares 
(zero  intercept)  and,  under  the  circumstances  of  deriva- 
tion, can  only  be  represented  as  a  reasonably  strong 
hypothesis  for  the  relation  at  hand.  Validation  and/or 
model  improvement  must  be  left  to  future  studies. 
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The  effects  of  precipitation,  growing  season,  and  soil 
potassium  on  forage  production,  as  expressed  by  this 
model,  are  shown  in  figure  1.  (See  equation  1,  appendix 
A.)  As  might  be  expected,  figure  1  shows  that  low  precipi- 
tation is  limiting  to  plant  growth,  regardless  of  length  of 
growing  season  (GS)  or  level  of  potassium  (K).  In  semiarid 
areas  of  Arizona  and  New  Mexico  where  precipitation 
(PR)  varies  between  5  and  10  inches  (13  and  25  cm)  and 
GS  is  generally  in  excess  of  120  days,  production  aver- 
aged from  only  100  to  250  lb/acre  (112  to  280  kg/ha)  on 
undisturbed  sites. 

With  PR  in  the  10-  to  15-inch  (25-  to  38-cm)  range  and 
GS  varying  from  perhaps  60  to  130  days,  typical  of  Wyom- 
ing and  Montana  mining  sites,  production  generally 
ranged  from  250  to  1,600  lb/acre  (280  to  1  793  kg/ha)  in 
the  presence  of  relatively  high  K  levels  (260  p/m).  Lesser 
amounts  of  K  (150  p/m)  were  associated  with  a  produc- 
tion drop  of  100  to  300  lb/acre  (112  to  336  kg/ha)  within 
the  60-  to  95-day  range  in  GS. 

North  Dakota  mining  sites  were  in  a  slightly  higher  PR 
zone  (16  to  17  inches  [41  to  43  cm])  and  had  moderately 
long  growing  seasons  (113  to  129  days).  Production  esti- 
mates for  this  region  of  the  model  were  heavily  affected 
by  the  more  copious  Wyoming  and  Montana  data  and 
only  reached  the  800-  to  900-lb/acre  (897-  to  1  009-kg/ha) 


range,  about  300  to  400  lb/acre  (336  to  448  kg/ha)  less 
than  the  average  of  the  few  observations  available  from 
the  undisturbed  North  Dakota  sites.  This  fact,  along  with 
the  relatively  high  productivity  rates  generally  accredited 
to  the  northern  Great  Plains,  suggests  the  need  for 
stronger  data  from  that  area.  Note  that  the  potassium 
effect  here  was  not  detectable  at  GS  >  100  days. 

A  few  higher  elevation  (7,000  to  8,000  ft  [2  134  to 
2  438  m])  mine  sites  with  precipitation  in  the  16-  to  23- 
inch  (41-  to  58-cm)  range  and  shorter  (50  to  81  days)  grow- 
ing seasons  were  sampled  in  northwestern  Colorado. 
Production  is  generally  greater  there,  ranging  from  about 
1,600  to  2,100  lb/acre  (1  793  to  2  354  kg/ha)  for  K  >  260 
p/m,  dropping  300  to  500  lb/acre  (336  to  560  kg/ha)  for  K 
<  150  p/m. 

As  a  whole,  this  model  provides  our  best  estimates  of 
minimum  production  requirements  or  standards  against 
which  the  success  of  planned  revegetation  efforts  on 
areas  to  be  mined  can  be  evaluated. 

Examples  of  estimates  of  forage  production  made 
from  the  equation  of  this  model  for  selected  values  of 
precipitation,  growing  season,  and  soil  potassium  are 
presented  in  table  1  to  demonstrate  the  output  format 
from  the  user-oriented  computer  program  for  this  model, 
which  is  shown  in  appendix  A. 


NUAL 
TATION 
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GROWING  SEASON  LENGTH 
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Figure  1. — Effects  of  annual  precipitation,  growing  season 
lengtfi,  and  soil  potassium  content  on  forage  pro- 
duction of  unmined  sites  on  surface  mineable  coal 
lands. 
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Table  1.— Effects  of  annual  precipitation,  growing  season  length,  and  soil  potassium  content  on  forage  production  of  unmined  sites 
on  surface  mineable  coal  lands  (lbs/acre) 
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Plant  Cover  Density  Model  for  Unmined  Areas 

Procedures  for  modeling  the  density  of  plant  cover  on 
unnnined  areas  were  similar  to  those  used  for  forage  pro- 
duction. Screening  of  linear  effects  again  revealed  that 
annual  precipitation  and  length  of  growing  season  were 
strong  variables  (Pr  <  0.025)  and  that  potassium  showed 
strength  at  short-to-medium  length  growing  seasons  and 
medium-to-high  precipitation  ranges.  An  interactive 
model  involving  these  variables  was  generated  and  fitted 
to  data  from  the  83  transects  located  on  unmined  areas. 
Relations  between  plant  cover  density  and  annual  pre- 
cipitation, growing  season  length,  and  soil  potassium 
content  are  illustrated  by  the  response  surfaces  in  figure 
2.  (See  equation  2,  appendix  A.)  As  expected,  this  model 
shows  that  low  PR  is  limiting  to  plant  cover  density 


development.  This  is  accentuated  somewhat  for  the  gen- 
erally longer  and  hotter  growing  seasons  of  the  semiarid 
Southwest  where,  with  5  to  10  inches  (13  to  25  cm)  of 
PR,  plant  cover  density  ranges  from  near  zero  to  25  per- 
cent. 

In  Wyoming  and  Montana,  at  10  to  15  inches  (25  to  38 
cm)  of  PR  and  60-  to  120-day  growing  seasons,  density  of 
cover  varies  from  30  to  85  percent.  This  full  range  appar- 
ently occurs  in  response  to  PR  alone,  where  growing 
seasons  are  in  excess  of  100  days.  In  northwestern  Colo- 
rado at  higher  elevations  (7,000  to  8,000  ft  [2  134  to 
2  438  m]),  with  fewer  GS  (50  to  81)  days  and  relatively 
high  K  levels  (  >  260  p/m),  the  density  of  plant  cover 
ranges  from  75  to  85  percent.  A  decrease  of  60  p/m  of 
potassium  occasions  a  reduction  of  about  10  percent  in 
cover  here. 


50  70  90  110  130  150 

GROWING  SEASON  LENGTH 
(Days) 


Figure  2. — Effects  of  annual  precipitation,  growing  season 
length,  and  soil  potassium  content  on  vegetative 
cover  density  of  unmined  sites  on  surface  mineable 
coal  lands. 
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Table  2.— Effects  of  annual  precipitation,  growing  season  length,  and  soil  potassium  content  on  vetetative  cover  density  of  unmined 
sites  on  surface  mineable  coal  lands  (percent) 
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At  a  PR  of  about  15  inches  (38  cm),  GS<110,  and 
K  >  260  p/m,  the  maximum  density  potential  of  plant 
cover  appears  to  be  about  85  percent,  with  a  drop  to  75 
percent  at  lower  K  levels  ( <  200  p/m).  For  GS  >  110 
days,  an  almost  full  cover  of  98  percent  is  reached  at 
PR  >  18  Inches  (46  cm).  These  conditions  are  typical  of 
the  North  Dakota  area.  Under  these  conditions,  potas- 
sium did  not  seem  to  be  a  limiting  factor. 

The  equation  for  this  model  was  used  to  estimate  plant 
cover  densities  for  selected  values  of  annual  precipita- 
tion, growing  season  length,  and  soil  potassium  content. 
Examples  of  these  estimated  plant  cover  densities  are 
shown  in  table  2;  they  demonstrate  the  output  format 
from  the  user-oriented  computer  program  for  this  model, 
which  is  shown  in  appendix  A. 

This  model  provides  estimates  of  plant  cover  density 
on  unmined  areas.  These  estimates  can  be  considered  as 
minimum  requirements  or  standards  against  which  the 
success  of  revegetation  efforts  on  mined  areas  can  be 
evaluated. 


Forage  Production  Models  for  Mined  Areas 

Two  models  were  developed  to  estimate  forage  pro- 
duction to  be  expected  from  revegetation  of  surface- 
mined  areas.  One  model  was  generated  to  estimate  the 
amounts  of  forage  produced  on  revegetated  areas  domi- 
nated by  native  plant  species.  Another  was  developed  to 
provide  similar  estimates  where  revegetated  areas  are 
dominated  by  introduced  plant  species.  Each  model  is 
composed  of  an  interactive  base  and  a  group  of  additive 
effects. 

The  age  of  the  planting,  the  climatic  components  of 
annual  precipitation  and  growing  season  length,  and  the 
potassium  content  of  spoils  (fitted  in  that  order)  added 
significantly  (Pr  <  0.10)  to  linear  regressions  for  forage 
production  from  both  native  and  introduced  vegetation. 
Subsequently,  attention  was  focused  on  interactive  rela- 
tions between  forage  production  and  the  three  indepen- 
dent variables:  age,  annual  precipitation,  and  growing 
season  length.  Forage  production  was  expected  to  be 
convex  upward  over  age,  reaching  a  maximum  at  some 
point  in  time;  to  be  concave  upward  over  precipitation; 
and  to  reach  a  peak  somewhere  within  the  broad  range  of 
growing  season  lengths  encountered.  Expected  trends 
were  adjusted  to  the  partitioned  data  graphically,  coordi- 
nated to  arrive  at  the  interactive  relation,  and,  simultane- 


ously, formulated  mathematically.  The  resulting  model 
was  adjusted  to  the  mined-transect  source  data  for 
native  (n  =  44)  and  introduced  (n  =  33)  species.  Under 
the  circumstances  of  derivation,  these  models  can  only 
be  represented  as  reasonably  strong  hypotheses  for  the 
relations  at  hand.  Model  validation  and/or  improvement 
must  be  left  to  future  studies.  The  interactive  portion  of 
this  model  for  forage  production  from  native  plant 
species  is  illustrated  by  the  response  surfaces  in  figure  3. 
(See  equation  3,  appendix  A.) 

Production  for  this,  the  interactive  portion  of  the 
model,  appeared  to  reach  an  upper  asymptote  at  age-of- 
planting  =  5  years.  Production  trends,  as  yet  unadjusted 
for  the  effects  of  potassium,  sodium,  pH,  and  vegetative 
rehabilitation  treatments,  still  closely  reflected  the  shape 
of  the  PR  effect  in  the  production  model  for  unmined 
areas  and  the  optimum  at  85  days  for  the  GS  effect  still 
existed. 

A  comparable  model  for  estimating  forage  production 
on  revegetated  areas  dominated  by  introduced  plant 
specis  is  identical  in  shape  to  the  model  for  native  plant 
species.  These  models  of  well-defined  interactive  effects 
were  refitted  to  their  respective  data  sets  and  were 
adopted  as  fixed  prediction  bases.  The  residuals  from 
these  bases  were  then  expressed  as  linear  additive 
effects  of  the  somewhat  weaker  continuous  variables, 
potassium,  sodium,  pH,  and  of  the  discrete  (present  or 
absent,  in  this  case)  revegetation  treatment  variables, 
tilling,  seeding  method,  topsoiling,  fertilizing,  supple- 
mental irrigation,  mulching,  and  seeding  time.  The  sum 
of  the  interactive  and  linear  additive  effects  constitutes 
the  estimate  of  production. 

Selected  values  of  annual  precipitation  and  growing 
season  length  were  utilized  in  the  equation  for  the  inter- 
active model  shown  in  figure  3.  The  output  from  this 
equation  together  with  the  additive  effects  of  selected 
values  of  soil  potassium,  sodium,  and  pH  content,  and 
those  of  the  seven  revegetation  treatments,  provided 
estimates  of  the  amounts  of  forage  produced  at  5  years 
of  age  on  surface-mined  areas  revegetated  with  predomi- 
nantly native  species.  Age  of  5  years  was  selected  for 
these  estimates  because  the  model  revealed  that  bio- 
mass  production  from  revegetation  of  western  coal  sur- 
face mines  tends  to  reach  peak  development  in  about  5 
years.  Examples  of  these  forage  production  estimates 
are  presented  in  table  3  to  demonstrate  the  output  format 
from  the  user-oriented  computer  program  for  this  model, 
which  is  shown  in  appendix  A. 
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Figure  3.— Effects  of  annual  precipitation,  growing  season  length, 
and  age  of  planting  on  forage  production  of  mined  and 
revegetated  sites  on  coal  surface  mines. 


Information  in  this  sample  portion  of  model  output  is 
divided  into  three  blocks,  each  of  which  represents  a  dif- 
ferent level  of  annual  precipitation.  In  order  from  left  to 
right,  these  levels  are  5  inches  (13  cm),  15  inches  (38  cm), 
and  25  inches  (64  cm).  All  remaining  variables  are  identi- 
cal for  each  of  the  three  blocks.  The  growing  season 
length  is  85  days.  The  soil  potassium  content  is  200  p/m. 
The  soil  sodium  content  (Na)  is  300  meq/liter  and  the  soil 
pH  is  8.  Each  of  the  three  blocks  of  the  table  consists  of 
eight  columns,  seven  of  which  are  occupied  by  data  from 
the  seven  revegetation  treatments.  The  last  column  in 
each  block  records  the  dry  weight  of  forage  produced  in 
the  presence  of  revegetation  treatments  checked  with  the 
letter  X.  Tillage  of  the  soil  (TIL),  consisting  of  ripping, 
disking,  or  plowing  prior  to  seeding,  is  shown  with  an  X. 
A  blank  indicates  that  the  soil  was  not  tilled.  Drilling  as  a 
seeding  method  (SM)  is  indicated  by  an  X,  whereas 
broadcast  seeding  is  represented  by  a  blank.  Addition  of 
topsoil  prior  to  seeding  (TPS)  is  indicated  by  an  X;  a 
blank  shows  that  no  topsoil  was  added.  An  X  indicates 
that  fertilizer  (PER)  was  added;  a  blank  indicates  that  it 
was  not.  Supplemental  irrigation  (IRR)  is  indicated  by  an 


X;  no  irrigation  is  represented  by  a  blank.  Mulching  (MUL) 
is  indicated  by  an  X,  and  no  mulching  is  a  blank.  An  X  for 
seeding  time  (ST)  indicates  fall  seeding,  whereas  a  blank 
represents  spring  seeding. 

Forage  production  values  such  as  those  for  native 
species  in  table  3  can  be  compared  with  the  appropriate 
forage  production  values  for  unmined  areas  such  as 
those  in  table  1.  Production  values  from  mined  and  re- 
vegetated  areas  after  5  years  of  growth  that  are  equal  to 
or  greater  than  the  production  values  from  similar  un- 
mined areas  denote  successful  revegetation  in  relation  to 
production  of  vegetation  on  the  undisturbed  areas  which 
is  assumed  to  reflect  ecological  potentials. 

Similarly,  estimated  amounts  of  forage  produced  were 
generated  from  the  interactive  and  additive  components 
of  the  model  for  forage  production  of  introduced  species. 
These  production  values  are  12.5  percent  higher  than 
those  of  native  species.  This  is  not  surprising  since  plant 
species  introduced  for  range  improvement  in  the  West 
generally  have  been  selected  for  their  exceptional  hardi- 
ness and  productivity  characteristics. 


6 


QDiiiJ  I  incoojcucoor^CTOTioivco.—  oocJ'^ocoool^"T^D(»^D(»cJ'TW(£co■^^covo^^ocM<r5>— 
DLlJD^  1  O'— '^c\ii^vc\]coir)<ocoLO'TioiD<— »(i)co<r)in'Tino(7)Oii)r^(i>r^C\JOTooa>oivinr^ 
ce:d_o  i  u>lD<T■lno^•u3"^OTO)lr)C^^D'T01eo(oo^!^cucolr5cocoorv.o's^(nivco<o^ 


X  X 


X  XX 


X  CO 


X 


X 


X 


X  X 


X 


XXX 


T3 


<  o 
z  o 

CO 


o 


a: 


Ll! 
li_ 

0-5 
CL 
h- 

s: 


X 


X 


X 


XX 


X 


X  X 


X 


X 


X 


X  X  XXX  X 

X  X  X  X  X  X 

xxxxx  xxxxx 


XXX 


xxxx 


X  X 


xxxx 


xxxxx 


> 

0} 


£3_  OJ 


 + 


X 


X  X  X  X  X  X 


XXXXXXXXXXXXXXX 


■D 

c 


c 
o 

o 

3 
T3 

o 


03 

ffl 


c 
o 


o 
> 


< ' 


o 
D 

o 

D 


LlI 


o 

I 

LU 

I 

Q_ 
CO 


LU 


LlI  I 

LU  q;  1 

C_  O  I 

<  1 


X  CO 

CL 


<  o 
z  o 


o 
<\1 


CO  LO 
O  CO 


on  m 


I 

-!-  —  + 

i  I-  i 

I  CO  i 
j  I 
+  + 
I 


I 


<r)<nocDco'TC^i'TOO'-coinr^r^<r)iiOCOi/7in>— crs.—  cO'-coi^<r)ii)Ocvicoa3coooiDcv'^ 

^5  Xi  ^  CO  CO  O  CT)  O  OJ  O      ^  O  ^  NT  CO  ^  !/5  C\i  O  (\J      CD      C\J  O  Oj  ^  CT) 

■^vcooj»ij^C\!coco\rcoiococor^^<x>i^<— tv.^c\jr-^cocD'-o-Tcoc£)i^^coci3^r^'^Cvir^ 


3 

4-  + 

I  K  I 

!  ! 

I  —  I 


a:  I 

lU  1 


! 


I  CO 

i  Q_ 

I  H- 
+ 
I 

I  CO  1 
I  I 
+  + 
i  _!  1 
I  —  1 
I  !—  I 
 + 


X 


X 


X  XX 


X  X 


X 


X 


X 


X 


X 


X 


X 


xxx: 


XX 


X  XX 


X  X 


X 


xxx 


xxx 


xxxx 


X 


xxxxx 


xxxxx 


xxxx 


xxxxx 


xxxxxx 


XXXXXXXXXXXXXXX 


■D 

c 

ca 

•-  o 

t:  CO 

S  "55 

o  ^ 

u_  — 

Q.  CO 

—  9J 

'o  o 

Q.  en 

w  Q. 

2  -B 

«  c 

o  £ 

E  ^ 

o  c 

o  .E 

w  E 

LU  Q. 

I 

CO 

_a) 


Q  CSriiJ 

o  LU  q; 
d;  CL  o 


<  I 

■i-  —  -r 
1  I—  I 
I    CO  1 


<  o 
z.  o 

CO 


o 
o 


CO  lO 

CD  CO 


a:  in 


I 


ojiDCO'OOKOcDCDC\j'^ooo<7>>— ^-in^-^-■^ocou^colno^o^coln'— ooc5<— i^O!^<7>CDcooincjo^ 

Kco^    ■^o     ocv    ioc\i    cvojeoTCO'srinoj    cvi^cdk    ■^CvTCTiCD'^cDKCO    coeoKco  ir>'^<r)coir)cO'— 

T      T--  T--      'TCOr^'^CMh-CO  CO'TOJr^      '^CMCOr^^CMh-CO  COTOJIV  •^CJCDCO^CDCO 


1  _l  i 

1  D  i 

I  s:  1 

-i-  + 

i  ir  1 

I  ir  '■ 

+ 
1 


LU  1 
Ll  1 


i  CO  i 

!   CL  i 

1   !—  1 

f  + 


I 


i   CO  i 


X 


X  XX 


X  xxx 


X 


X 


X 


X 


X  X 


xxx  X 


xxx  X 


X  X 


X 


X 


xxx 


X 


xxxx 


X  xxx 


xxxx  xxxx 


X 


xxxxx 


xxxxx 


xxxxx 


xxxxxx 


XXXXXXXXXXXXXXX 


7 


Plant  Cover  Density  Models  for  Mined  Areas 

Procedures  for  modeling  the  density  of  plant  cover  on 
mined  areas  were  similar  to  those  used  to  model  forage 
production.  Two  interactive  models  were  generated,  one 
for  revegetated  areas  predominantly  characterized  by 
native  species  and  the  other  for  areas  dominated  by 
introduced  species.  The  models  relate  plant  cover  density 
to  the  amount  of  annual  precipitation,  age  of  the  plant- 
ing, and  length  of  the  growing  season.  These  relations 
for  mined  areas  revegetated  with  native  species  are  illus- 
trated by  the  response  surfaces  shown  in  figure  4.  (See 
equation  4,  appendix  A.) 

Plant  cover  density  increased  with  increasing  amounts 
of  annual  precipitation  and  with  increasing  age  of  vegeta- 
tion up  to  about  5  years.  It  was  also  greatest  at  about  85 
days,  an  intermediate  growing  season  length.  Selected 
values  of  annual  precipitation  and  of  growing  season 
length  were  used  in  the  equation  for  this  interactive 
model.  These,  along  with  the  additive  effects  of  selected 
values  of  soil  potassium,  sodium,  and  pH  content  and  of 
the  revegetation  treatments,  permitted  estimates  of  the 
plant  cover  densities  expected  at  5  years  of  age  on 
surface-mined  areas  revegetated  with  predominantly 
native  species. 


Examples  of  these  plant  cover  density  estimates  are 
presented  in  table  4  to  demonstrate  the  output  format 
from  the  user-oriented  computer  program  for  the  model, 
which  is  shown  in  appendix  A. 

Plant  cover  density  values,  such  as  those  in  table  4, 
can  be  compared  with  the  appropriate  cover  density 
values  for  unmined  areas,  such  as  those  in  table  2.  Plant 
cover  density  values  (age  =  5  years)  from  mined  and  re- 
vegetated areas  that  are  equal  to  or  greater  than  those 
from  similar  unmined  areas  denote  successful  revegeta- 
tion. Plant  cover  density  development  on  these  unmined 
areas  is  assumed  to  reflect  ecological  potentials. 

Similarly,  estimated  plant  cover  density  values  were 
generated  from  the  interactive  and  additive  components 
of  the  plant  cover  density  model  for  introduced  species. 
These  density  values  are  8.8  percent  lower  than  those  of 
native  species.  While  introduced  species  are  superior  in 
producing  forage  native  species  appear  able  to  provide 
better  protective  ground  cover.  Thus,  both  kinds  of  plants 
possess  characteristics  that  enable  them  to  meet  re- 
quirements of  the  Surface  Mining  Control  and  Recla- 
mation Act  by  "establishing  a  diverse,  effective,  and  per- 
manent vegetative  cover"  on  revegetated  coal  mine 
spoils  in  the  Western  United  States. 


AGE 


ANNUAL 
PRECIPITATION 
(Inches) 


150 


GROWING  SEASON  LENGTH 
(Days) 


Figure  4. — Effects  of  annual  precipitation,  growing  season 
length,  and  age  of  planting  on  vegetative  cover 
density  of  mined  and  revegetated  sites  on  coal 
surface  mines. 
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APPLICATION 

The  outputs  from  this  research  include  two  regression 
models  for  unmined  (undisturbed)  areas,  one  for  use  in 
estimating  forage  production  weights  and  the  other  for 
estimating  vegetative  cover  densities.  Also  included  are 
four  regression  models  for  mined  areas.  Two  models  are 
for  use  in  estimating  forage  production  weights,  one 
where  native  plants  predominate,  the  other  where  intro- 
duced plant  species  are  dominant.  The  other  two  models 
are  for  use  in  estimating  vegetative  cover  densities  of 
native  and  of  introduced  vegetation. 

Each  of  these  six  models  is  expressed  as  an 
algebraic  equation  in  appendix  A.  Following  each  of 
these  equations  are  FORTRAN  IV  computer  programs  to 
facilitate  computer  solutions  for  pertinent  combinations 
of  the  different  levels  of  the  independent  variables  in- 
volved and,  as  appropriate,  for  native  and  introduced 
species.  The  values  contained  in  tables  1,  2,  3,  and  4 
are  sample  outputs  from  these  FORTRAN  programs. 

Another  output  from  this  research  consists  of  two 
multicolored  maps  (scale  1:1,000,000)  that  depict  the  fol- 
lowing information  for  the  surface  mineable  coal  areas  of 
the  West,  as  well  as  for  most  of  the  intervening  areas 
(appendix  B,  inside  back  cover). 

1.  Annual  precipitation  (inches); 

2.  Growing  season  length  (days); 

3.  Natural  vegetation  types  and  their  identification 
numbers  (Kuchler  1964); 

4.  Soil  associations  and  their  identification  numbers 
(Aandahl  1972;  Cipra  and  others  1977;  Jay  and  others 
1975;  Maker  and  others  1974;  Southard  1973;  Wilson  and 
others  1975;  and  Young  and  Singleton  1977). 

These  maps  provide  estimates  of  precipitation  and 
growing  season  length  and  identify  the  soil  association 
and  natural  vegetation  type  for  any  location  on  surface- 
mineable  coal  areas  of  the  West. 

Here  are  instructions  covering  use  of  the  maps, 
equations,  and  tables  in  order  to  estimate  the  average 
degree  of  revegetation  success  on  any  of  these  surface- 
mineable  coal  areas. 

1.  For  any  selected  location  on  either  of  the  maps  de- 
termine values  of  precipitation  and  growing  season 
length  and  identify  the  soil  association  and  vegetation 
type  by  numbers.  The  soil  association  numbers  are  those 
contained  either  in  Soils  of  the  Great  Plains  (Aandahl 
1972)  or  in  individual  State  soil  description  publications 
(Cipra  and  others  1977;  Jay  and  others  1975;  Maker  and 
others  1974;  Southard  1973;  Wilson  and  others  1975;  and 
Young  and  Singleton  1977),  which  are  referenced  in  the 
Publications  Cited  section.  These  publications  contain 
information  concerning  the  chemical  and  physical  char- 
acteristics of  the  soil  associations.  The  vegetation  type 
numbers  are  contained  in  Kuchler's  1964  map  of  "The 
Potential  Natural  Vegetation  of  the  Conterminous  United 
States."  This  map  identifies  the  dominant  native  species 
that  make  up  each  vegetation  type.  Such  identification 
provides  ecological  guides  to  selection  of  plant  species 
adapted  for  revegetation  purposes. 

2.  Use  the  values  of  annual  precipitation,  growing 
season  length,  and  soil  potassium  content  in  equation  1, 
appendix  A  for  unmined  areas  to  calculate  forage  pro- 


duction in  pounds  per  acre.  Similarly,  the  values  of  these 
factors  can  be  used  in  equation  2,  appendix  A  to  calcu- 
late cover  density  in  percent.  These  values  are  con- 
sidered to  define  the  ecological  potential  for  producing 
forage  and  developing  cover  density  at  this  location.  If 
published  information  concerning  the  chemical  charac- 
teristics of  the  particular  soil  association  under  consider- 
ation is  not  available,  then  soil  sampling  may  be  neces- 
sary to  obtain  site  specific  data  on  the  soil  potassium 
content. 

3.  Depending  upon  whether  revegetation  of  mined 
areas  or  planned  revegetation  of  areas  to  be  mined  is 

with  native  or  with  introduced  types  of  vegetation,  enter 
the  appropriate  version  of  the  mined  areas  interactive 
component  equation  (equation  3,  appendix  A)  with  the 
same  precipitation  and  growing  season  values  that  were 
used  in  the  unmined  area  equations.  Calculate  the  inter- 
active effects  of  these  factors  on  the  estimated  pounds 
per  acre  of  forage  produced  at  age  5  years. 

4.  Enter  the  appropriate  additive  component  table  for 
forage  production  (table  1,  appendix  A).  Calculate  addi- 
tional effects  on  forage  production  of  spoil  potassium, 
sodium,  and  pH  conditions  encountered  or  expected,  as 
well  as  additional  effects  of  the  seven  revegetation  treat- 
ments indicated.  In  the  event  that  previous  mine  spoil 
analyses  are  not  available  or  do  not  provide  sufficient 
information  on  potassium,  sodium,  and  pH  levels,  spoil 
sampling  may  be  necessary  to  obtain  specific  on-site 
values  of  these  factors. 

5.  Add  the  pounds  per  acre  of  forage  production  calcu- 
lated from  the  interactive  component  equation  to  the 
pounds  per  acre  calculated  from  the  additive  component 
table.  These  values  are  the  estimated  forage  production 
amounts  achievable  under  the  climatic,  spoil,  and  revege- 
tation treatment  conditions  encountered. 

6.  Similar  calculations  of  achievable-  plant  cover  densi- 
ties can  be  made  utilizing  the  climatic,  spoil,  and  revege- 
tation treatment  factors  in  the  interactive  component 
equation  (equation  4,  appendix  A)  and  additive  compon- 
ent table  (table  2,  appendix  A)  for  plant  cover  density. 

7.  Select  those  alternative  revegetation  treatment  com- 
binations that  provide  forage  production  weight  and  plant 
cover  density  values  equal  to  or  greater  than  the  values 
from  the  comparable  unmined  area  equations.  These 
values  denote  success  compared  to  the  ecological 
standards  for  existing  conditions  and  they  identify  the 
combinations  of  climatic,  spoil,  and  revegetation  treat- 
ment conditions  needed  to  achieve  such  success  5  years 
after  planting. 

Following  is  an  example  of  the  use  of  the  maps, 
tables,  and  equations  to  compare  estimated  forage  pro- 
duction weights  of  predominantly  native  vegetation  on  a 
selected  unmined  and  mined  site  in  eastern  Montana. 

1.  From  the  maps,  it  is  determined  that  the  average 
annual  precipitation  at  the  selected  site  is  15  inches 
(38  cm).  The  growing  season  length  is  85  days.  The  soil 
association  is  Aandahl's  number  117,  which  is  moder- 
ately low  in  potassium  (less  than  200  p/m),  high  in 
sodium  (more  than  300  meq/L),  and  high  in  pH  (about  8). 
The  vegetation  type  is  Kuchler's  number  6,  grassland- 
sagebrush,  which  contains  a  number  of  highly  adapted 
and  suitable  grass  and  shrub  species  for  revegetation. 
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2.  Utilizing  this  information,  enter  table  1  in  the  center 
section:  200  p/m  of  potassium,  15  inches  (38  cm)  of  pre- 
cipitation, and  85  days  of  growing  season.  This  combina- 
tion of  conditions  results  in  a  tabulated  estimate  of  1,531 
pounds  per  acre  (1  716  kg/ha)  of  native  vegetation  weight 
produced  on  the  unmined  site,  the  ecological  potential 
for  this  site.  This  same  value,  1,531  pounds  per  acre  can 
be  obtained  by  solving  equation  1,  appendix  A  as  follows: 

N  =  1.8  -  0.56  *  e- 0-01770 
=  1.8  —  0.56  *  0.9825  =  1.250 

YPFL  =  1,570  +  1,060  *  e-O-OOO^^ 

=  1,570  +  1,060  *  0.9994  =  2,629.4 

YP1  =  YPFL  +  560  *  e  -  0-25025 

=  2,629.4  +  560  *  0.7786  =  3,065.4 

PRODUCTION  =  ^'065;4    ,      ^  25  *  0.94584 
(25)1.25 

=  1,531  pounds/acre 

3.  Further  utilizing  the  information  obtained  from  the 
maps  and  the  cited  soil  reference  (Aandahl  1972),  enter 
table  3  in  the  center  section.  Here,  the  native  vegetation 
production  weight  estimates  are  for  15  inches  (38  cm)  of 
precipitation,  85  days  of  growing  season,  200  p/m  of 
potassium,  300  meq/liter  of  sodium,  and  pH  of  8— the 
conditions  encountered  at  this  site.  This  combination  of 
conditions,  together  with  the  various  revegetation  treat- 
ment combinations  indicated  in  the  table,  results  in  a 
number  of  production  weight  estimates  in  excess  of 
1,531  pounds  per  acre  (1  716  kg/ha)  when  the  vegetation 
on  the  mined  areas  is  5  years  old.  The  highest  production 
weight  in  that  portion  of  the  total  table  represented  by 
table  3—2,214  lb/acre  (2  482  kg/ha)— can  be  achieved  by 
selecting  a  combination  of  treatments  that  includes  drill- 
ing seed  in  the  spring,  fertilizing  the  area  with  a  well- 
balanced  N— P— K  fertilizer,  and  irrigating.  Very  nearly 
the  same  production  level  (2,186  lb/acre  or  2  450  kg/ha) 
can  be  achieved  by  tilling,  drilling  in  the  fall,  and  fertil- 
izing without  irrigation.  Almost  as  much  production  (2,173 
lb/acre  or  2  436  kg/ha)  is  obtainable  with  much  less  effort 
and  expense  simply  by  drilling  and  fertilizing  in  the  fall 
without  either  tillage  or  irrigation. 

These  estimates  of  forage  production  weight  can  be 
obtained  by  solving  equation  3,  appendix  A  and  adding 
the  appropriate  elements  from  the  additive  component 
table  (table  1,  appendix  A).  This  solution  for  the 
maximum  weight  of  2,214  lb/acre  (2  482  kg/ha)  is  as 
follows: 
a.  Equation  3 

PRODUCTION  =  0.0061896  *  YPPR  *  YPGS 

*  (PR) "I -6  *  1.04368 

YPPR  =  e- 0.00510  ^  0.99491 

YPGS  =  940  +  2,510  *  e-0 

=  940  +  2,510  *  1  =  3,450 

PRODUCTION  =  0.0061896  *  0.99491  *  3,450 

*  (15)'' -6  *  1.04368 
=  1,688  pounds/acre 


b.  Additive  components  (native  vegetation) 

Soil  potassium  (K)  200  *  5.4         =  1,080 

Soil  sodium  (Na)  300  *  (-0.088)  =  -26 

Soil  acidity  (pH)  8*117.9      =  944 

Intercept  -2,216 

Seed  drilled  341 

Fertilizer  369 

Irrigation  35 

Spring  seeding  0 


Additive  components  total         526  pounds/acre 
Total  production  (1,688  -i-  526)     2,214  pounds/acre 

c.  These  forage  production  weight  values— up  to  2,214 
pounds  per  acre  (2  482  kg/ha)  on  revegetated  mined  areas 
compared  to  only  1,531  pounds  per  acre  (1  716  kg/ha)  on 
similar  but  unmined  sites— indicate  that  successful  re- 
vegetation  can  be  achieved  on  mined  sites  characterized 
by  the  climatic,  soil,  and  treatment  conditions  specified 
in  this  example.  Still  unanswered  are  questions  concern- 
ing the  permanency  of  production  weights  of  such  mag- 
nitude in  the  absence  of  further  fertilization  or  of  what 
the  land  management  requirements  are  to  maintain  such 
production.  Examination  of  the  complete  computer- 
produced  tables  leads  one  to  conclude,  however,  that  re- 
vegetation  treatment  alternatives  are  available  to  pro- 
duce forage  weights  (and  cover  densities)  which,  5  years 
after  establishment,  are  greater  than  the  ecological 
potentials  of  most  western  coal  mine  sites. 

4.  The  algebraic  specifications  for  the  six  models  are 
presented  in  appendix  A  as  equations.  Immediately  fol- 
lowing each  is  an  associated  FORTRAN  IV  computer 
program  designed  to  produce  tables  like  the  text  examples 
(tables  1-4).  Note  that  these  programs  are  designed  to  run 
on  IBM  or  IBM-compatible  systems  such  as  the  Amdahl 
470V/6-II.  Minor  program  changes,  such  as  "read"  and 
"write"  instructions  and  a  substitute  for  the  standard 
IBM  ERRSET  subroutine,  may  have  to  be  made  before 
they  will  run  on  alternative  systems. 

In  the  programs  for  mined  areas,  three  groups  of  pre- 
selected values  of  PR,  GS,  K,  Na,  and  pH  (in  this  se- 
quence) are  used  in  statement  003  as  a  major  control  for 
program  output.  The  values  are:  5,  50,  0,  0,  4;  15,  85,  200, 
150,  6;  and  25,  120,  400,  300,  8.  These  can  be  changed  to 
accommodate  alternative  output  needs,  but  must  be 
within  the  limits  of  use  specified  for  each  variable  and 
algebraic  model.  The  output  format  for  vegetative  treat- 
ments is  considered  fixed  and  should  not  be  altered. 

The  program  for  unmined  areas  produces  both  forage 
production  and  cover  percent  estimates.  Only  three  vari- 
ables control  the  output  here,  PR,  GS,  and  K.  Three  pre- 
selected values  of  each  appear  in  the  program  as 
follows: 

Statement  0008,  where  K  =  50,  120,  35.  K  represents 
GS  here  and  the  statement  is  interpreted  to  mean  that 
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GS  levels  start  at  50,  and  are  Incremented  by  35  until  the 
maximum  of  120  is  reached.  Thus,  GS— levels  =  50,  85, 
and  120. 

Statement  0010,  where  LL  =  1,  401,  200.  LL  represents 
K  and  the  statement  is  interpreted  as  In  0008,  so  that  K 
levels  are  1,  201,  and  401.  "One"  has  been  added  to  the 
desired  values  of  0,  200,  and  400  to  accommodate  a  pro- 
gramming limitation  at  zero.  One  is  subtracted  (in  state- 
ment 0011)  before  model  computations  are  begun  so  that 
the  final  output  is  for  K  levels  of  0,  200,  and  400. 

Statement  0012,  where  M  =  5,  25,  10.  M  represents  PR 
and  the  statement  is  interpreted  as  in  0008,  so  that  the 
PR  levels  are  5,  15,  and  25. 

The  levels  of  these  variables  can  be  changed  to  meet 
alternative  output  needs  but,  again,  must  be  within  the 
limits  specified  for  each  variable  and  algebraic  model. 

Neither  pagination  nor  indexing  is  provided  in  these 
programs.  It  is  suggested  that  users  add  such  provisions 
for  convenience  in  application  of  the  many  tables  gener- 
ated (243  tables  here). 


DISCUSSION 

The  primary  objectives  of  this  investigation  were  to 
develop  capabilities  for  estimating  the  degree  of  revege- 
tation  success  to  be  expected  under  a  wide  variety  of 
climatic  conditions,  soil  and  spoil  properties,  and  revege- 
tation  treatments  utilizing  site-specific  revegetation  data 
and  information  from  most  of  the  coal  surface  mines  in 
the  interior  West.  The  more  important  study  develop- 
ments are  these: 

•  A  strong  conceptual  framework  for  evaluating  the 
success  of  proposed  vegetative  rehabilitation  efforts 
on  areas  to  be  surface  mined. 

•  Site-specific  maps  have  been  developed  to  provide 
precipitation,  growing  season  length,  soil,  and  vege- 
tation type  information,  critical  to  the  evaluation 
system. 

•  Results  of  most  existing  surface-mine  rehabilitation 
efforts  in  the  interior  West  through  1976  were  sur- 
veyed as  a  basis  for  evaluating  revegetation 
success. 

•  Interim  predictors  (models)  were  developed,  based 
on  these  survey  data,  for  forage  production  and 
cover  density  potentials  on: 

•  Undisturbed  sites  adjacent  to  the  study  mines. 

•Mined  areas  with  various  combinations  of  planting 
and  postplanting  methods  and  treatments  designed 
to  enhance  the  total  rehabilitation  effort  for: 

•  Native  species; 

•  Introduced  species. 

This  investigation  shows  that  practical  criteria  for 
measuring  revegetation  success,  namely,  the  amount  of 
forage  produced  and  the  density  of  plant  cover  devel- 
oped, are  affected  significantly  by  at  least  two  major  cli- 


matic factors  that  are  not  readily  susceptible  to  alter- 
ation (PR  and  GS);  by  three  properties  of  spoil  materials 
that  are  subject  to  limited  modification  through  manage- 
ment (K,  Na,  and  pH).  Revegetation  success  is  also 
influenced  by  seven  revegetation  treatments,  each  of 
which  provides  at  least  two  management  alternatives, 
and  by  the  age  of  the  vegetation.  These  characteristics 
account  for  about  one-half  to  three-fourths  of  the  total 
variance  in  forage  production  and  plant  cover  density  in 
the  prediction  models  (table  3,  appendix  A). 

It  is  emphasized  here  that  the  models  should  be  used 
with  considerable  discretion  since  the  source  data,  while 
the  best  available,  is  weak  in  many  respects.  First,  there 
were  only  34  mines  in  the  entire  area  covered  by  this 
study  where  revegetation  was  even  attempted.  The  kinds 
of  revegetation  efforts,  insofar  as  specific  treatments  or 
combinations  of  treatments  are  concerned,  differ 
between  mines  and,  in  many  cases,  even  between  years 
on  the  same  mine.  For  example,  seeds  of  native  species 
may  be  broadcast  in  the  fall  on  the  spoils  of  one  mine, 
whereas  seeds  of  introduced  species  may  be  drilled  into 
the  soils  of  another  mine  in  the  spring.  In  each  case,  the 
combinations  of  other  treatments,  such  as  topsoiling, 
tilling,  mulching,  or  fertilizing,  might  be  substantially 
different.  Of  the  large  number  of  possible  treatment 
combinations,  relatively  few  exist  on  the  28  mines 
studied,  and  the  comparative  revegetation  success  of 
those  that  do  exist  is  confounded  by  mine-to-mine  differ- 
ences in  climatic  and  growing  media  environment. 
Further,  "treatments"  cannot  be  considered  to  be  stan- 
dardized. For  example,  the  depth  and  quality  of  topsoil 
(when  added)  are  almost  certain  to  differ  between  mines. 
Fertilizer  composition  and  rates  and  times  of  application 
are  largely  unknown.  Also  unknown  are  rates  of  seed 
application,  depths  of  drilled  seed,  the  quality  of  seed, 
the  mix  of  species  seeded,  the  depth  of  tillage,  the  kinds 
and  amounts  of  mulch  applied,  and  supplemental  irriga- 
tion amounts  and  frequencies.  At  the  time  of  this  study, 
about  half  of  the  revegetation  efforts  were  no  more  than 
2  years  old;  so,  on  the  mines  involved,  seeded  plants  had 
not  had  time  to  respond  fully  to  the  growing  environment 
and  revegetation  treatments.  Revegetation  success,  of 
course,  is  strongly  affected  by  the  particular  amount  and 
timing  of  precipitation  during  the  year  or  period  of  years 
involved  in  the  revegetation  effort. 

It  seems  certain  that  the  interest  in  western  coal  de- 
velopment and  associated  environmental  controls  will 
continue  for  a  number  of  years.  Under  these  circum- 
stances and  because  of  the  low-order  quantitative  infor- 
mation available  on  the  success  of  spoil  rehabilitation 
efforts,  we  recommend  strongly  that  extensive,  well- 
planned  experimentation  be  underwritten  by  the  agencies 
requiring  this  information.  Although  such  research  is  now 
underway  in  the  Forest  Service,  it  is  of  limited  scope  and 
cannot  fully  meet  all  of  the  management  needs  specified. 
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APPENDIX  A—  EQUATIONS  AND 
SUPPLEMENTAL  TABLES  FOR  USE  IN 
ESTIMATING  FORAGE  PRODUCTION  WEIGHTS 
AND  PLANT  COVER  DENSITIES 

The  figures  and  partial  tables  presented  in  this  paper 
are  useful  as  an  aid  in  understanding  the  relations  de- 
scribed by  the  models  for  estimating  forage  production 
weights  and  plant  cover  densities.  They  are  not  particu- 
larly effective,  however,  as  tools  for  generating  such 
estimates.  Needed  to  estimate  useful  values  of  forage 
production  and  plant  cover  density  are  the  equations  for 
the  interactive  portions  of  these  models  (both  unmined 
and  mined  area  models)  and  the  tables  for  the  additive 
components  (mined-area  models  only).  These  equations 
and  additive-effect  tables  follow.  Note  that  the  sum  of 
the  interactive  fixed  base  and  linear  residual  effects 
constitutes  the  estimate  of  total  forage  production 
weight  or  vegetation  cover  density.  While  these  models 
can  be  solved  readily  with  most  pocket  calculators,  com- 
puterized solutions  are  strongly  recommended  if  produc- 
tion of  tables  of  estimates  is  the  desired  end  product. 


1.  Equations  for  unmined  areas 

The  models  for  unmined  areas  are  basically  inter- 
active. 

a.  Forage  production  model  equation       (Equation  1) 


Production  = 


(y^(pr)n\ 


N  =  1.8  -  0.56  *  e 
If  GS  <  87 


*  0.94584 


(180  -  GS) 
106 


-  1 


0.3 


YP1  =  YPFL  +  560  *  e 


K 

450 


0.6 


GS 
86.6 


-  1 


YPFL  =  1570  +  1060  *  e 


0.12 


If  GS  >  87 


YP2  =  1700  +  (YPI^-  1700)  *  e 


(180-GS) 
93.4 
0.09 


LIMITS 

5  <  PR  <  25  PR  =  annual  precipitation,  inches 

50  <  GS  <  180  GS  =  growing  season,  days 

0  <    K  <  450  K  =  potassium,  p/m 


@  YPFL  =  2630 
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b.  Plant  cover  density  model  equation  (Equation  2) 

PR     ^  ' 


Cover  =  YP  * 


25 


(1-1) 


1  -e         -  I)} 


*  0.97917 


YP  =  YPFL  +  YPAD 


GS 
180 


-  1 


YPFL  =  80  +  20  *e 


0.46 


K 

400 


-  1 


GS 
180 


-  1 


YPAD  = 


9  *  e 


0.43 


1  -  e 


0.46 


GS 
180 


-  1 


I  =  0.14  +  0.285  *  e 


0.46 


LIMITS 

5  <  PR  <  25  PR  =  annual  precipitation,  inches 

50  <  GS  <  180  GS  =  growing  season,  days 

0  <    K  <  450  K  =  potassium,  p/m 
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c.  FORTRAN  IV  computer  programs  for  estimating  forage  production  (equation  1) 
and  plant  cover  density  (equation  2) 


0001 
0002 
0003 
0004 
0005 
0006 
0007 
0008 
0009 
001  0 
001  1 


001 
001 
001 
001 
001 
0017 
001  8 
001  9 
0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 
0028 
0029 
0030 
0031 
0032 
0033 
0034 
0035 
0036 
0037 
0038 
0039 
0040 
0041 
0042 


0043 
0044 


0045 
0046 
0047 


1  8 
20 

25 


26 

27 
28 
30 
100 

101 
102 


REAL  I,LN 

DIMENSiaN  CLP (9) 

CALL  ERRSET  (208,256,-1,0) 

PRINT  101 

DS  30  J=1 , 2 

IF(J,EQ.2)   PRINT  100 

PRINT  102 

DO  25  K=50, 1 20, 35 

N  =  0 

DO  20  LL=1 , 401 , 200 
L=LL-1 

DS  18  M=5, 25, 10 

GS  =  K 

PR=M 

PQ  =  L 

N  =  N+1 

IF( J, EQ. 2)  GQ  TQ  26 

I  =  , 1 4+ , 285* (EXP( -(ABS( (GS/1 80 , - 1  , ) / , 46 ) « *  1 5 ) ) ) 
YPFL=80, +20, * (EXP( "(ABS( (GS/1 80 , - 1  , ) / , 46 ) *  *  1 5 ) ) ) 
YPA=1 , -(EXP(-(ABS( (GS/1 80, -1,)/,45)««15))) 
YPD=9, *(EXP(-(ABS( (PQ/400. -1,)/,43)«*15))) 
YP=YPFL+YPA*YPD 

LN=EXP( -(ABS( (PR/25. -1 .)/(1,-I))*^5n 
RN  =  EXP(  -(  ( 1  ,  /( 1  ,  -n  )  «*5)  ) 
CLP(N) = ( (LN-RN)/( 1 , -RN) ) *YP* 
CQNTI NUE 
CGNTI NUE 

PRINT  103, K, (CLP(II), 11=1,3) 
CQNTI NUE 
PRINT  104 
GG  TQ  30 

YN=1  , 8-,56*(EXP(-(ABS( ( ( 1 80 , -GS )/1 06 , -1  ,  )/ , 3) **3 , 

YPFL=1570, +1060. *(EXP(-(ABS( (GS/86, 6-1 . )/, 12)**4) 

YP1 =YPFL+560. *(EXP(-(ABS( (PG/450, -1,)/,6)**18))) 

I F(GS , GT, 86, 5)   GG  TG  28 

CLP(N)  =  (  (YP1/25  ,  **YN)  *PR**YN)  =^  ,  94584 

GG  TQ  18 

+  (YP1-1700.  )*(EXP(-(ABS( (  (180 


97917 


K,  (CLP( I  I ) ,  I  I =4, 6) ,K,  (CLPCI I ) ,  I  I =7, 


8)  )  ) 

)  ) 


YP1=1700 
GQ  TG  27 
CQNTI NUE 
STQP 

FQRMAT( ' 1 ' 
*'F  G  R  A  G 
*'  PER 

FORMAT ( ' -' 
*'F  0  R  A  G 

FGRI1AT(  '0' 


-GS)/93,4-1 , )/,09)*«3))) 


39X, 'U 
EV 

A  C  R 
39X, 'U 
EV 
7X, 'NO 


N  n 
28X, 
E'  ) 
N  M 
34X, 


NED 
PROD 

NED 
P  E  R  C 


ARE 

U  C  T  1 

ARE 
E  N  T  / 


A  SV 
Q  N 

A  SV 
,  G  E 


0 


47X, 
0  u 

47X, 
F 


N  D  S 


COVER') 


POTASSIUM' , 20X, '200  PARTS/M I LL I GN  POTASSIUM' 


*12X, '400  PARTS/MILLION 
*3( '+', 7( '='),'+', 17( ' - 
«3( ' ; DAYS  OF 
*3( ' : GROWING 
*3( ' ;   SEASON;     5  1 5 

*3( '+', 7( '-'),'+', 1 7( ' -' 

103  FORMAT( '0' , 3( ' ;  ',13,' 

104  FGRMAT( ' 0' , 3( ' +' , 7( ' -' ) 
END 


POTASSIUM'/' 
),'+', 12X)/' 


INCHES 
PRECIPITATION 
25 
)  '+' 


, 12X)/' 
, 12X)/' 
, 12X)/' 
, 1 2X) ) 
2(F5, 0, IX) 
17( ' -' )    ' + 


F5, 0, ' 
, 12X) ) 


12X)  ) 
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2.  Equations  for  mined  areas 

The  nnodels  for  mined  areas  are  partially  interactive 

and  partially  additive. 

a.  Forage  production  model 

(1)  Interactive  component  equation       (Equation  3) 

Kind  of 
vegetation 

ProCi^on  =  0.0061896  *  YPPR  *  YPGS  *  (PR)^  ^  *  i?7448  (Soduced) 


YPPR  =  e 


AGE 
7 

0.9 


If  GS  <  85 


YPGS  =  940  +  2510  *  e 


GS 

85  -  ^ 
0.16 


If  GS  >  85 


GS 

85  

0.12 


-  1 


YPGS  =  2250  +  1200  *  e 
LIMITS 


5  <  PR    <  25      PR  =  annual  precipitation,  inches 
50  <  GS    <  180     GS  =  growing  season,  days 
0  <  AGE  <  7      AGE  =  age  of  planting,  years 


(2)  Additive  component  table  (Table  1) 

Additive  forage  production 
Component   (lbs/acre)  


Native 

Introduced 

Soil  potassium  (K) 

5.4 

'+ 

4.1 

Soil  sodium  (Na) 

0.088 

1  

1.4 

Soil  pH 

117.9 

143.4 

Tillage  (Til) 

+ 

13.0 

308.0 

Seeding  method  (SM) 

broadcast 

0 

0 

drilled 

+ 

341.0 

783.0 

Topsoil  (TPS) 

83.0 

+ 

99.0 

Fertilizer  (PER) 

+ 

369.0 

+ 

430.0 

Irrigation  (IRR) 

+ 

35.0 

+ 

418.0 

Mulching  (MUL) 

178.0 

+ 

160.0 

Seeding  time  (ST) 

spring 

0 

0 

fall 

+ 

335.0 

78.0 

Intercept 

3_ 

2216.0 

3  _ 

1181.0 

'Multiply  these  production  values  by  parts  per  million  of  soil  component 
(potassium  or  sodium). 

'Multiply  these  production  values  by  pH  units  of  soil  acidity. 
^These  intercept  values  must  be  added  to  the  total  value  of  forage  pro- 
duction. 

LIMITS 


0  <  K    <  450 
0  <  Na  <  1000 
4  <  pH  <  9 
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(3)  FORTRAN  IV  computer  programs  for  estimating  forage  production  (equation  3)  for: 
(a)  Native  vegetation 


0001 

DIMENSION  IX(5, 3) , LB(3) , A(7)  ,XT(3) 

0002 

INTEGER*2  N(7) 

0003 

DATA  IX/5, 50, 0, 0, 4, 15, 85, 200, 150, 6, 25, 120, 400, 300, 8/, 

*N/7*0/,X/'X'/, A/7*'    '/,B/'  '/ 

0004 

CALL  ERRSET  (208,256,-1,0) 

0005 

DO  45  Mi^=l  ,  3 

0006 

DO  45  M=1 , 3 

0007 

DO  45  L=1 , 3 

0008 

DO  45  K=l,3 

0009 

DO  20  1=1,3 

001  0 

Z=IX(2,K) 

001  1 

YPPR= (EXP( -(ABS( (5, 0/7 ,0-1 .0)/.9)**4.6))) 

0012 

IF( IX(2,K) .GT. 85)  GO  TO  1 0 

0013 

YPG1 =2510«(EXP(-(ABS( (Z/85. -1,)/.16)**4)) )+940 

0014 

XT( I ) =YPPR^YPG1 * , 0061 8959* FLOAT ( IX(1,I))**1.6 

001  5 

GO  TO  15 

001  6 

10 

YPG2=1  200* (EXP( -(ABS( (Z/85, -1,)/,12)**4)) )+2250 

001  7 

XT (I ) =YPPR«YPG2* , 0061 8959xFL0AT( IX( 1 , I ))**1 .6 

001  8 

15 

XT( I )=XT( I )«1 , 04368 

c  *** 

FOR  ZERO  COMPUTATION,  *** 

0019 

LB( I )=-2215. 925+XT( I ) +IX(3, L ) *5 , 39904+ I X( 4 , M) * ( - , 08788242) + 

*IX(5,MM)*1 17. 9473 

0020 

IF(LB( I ) . LT. 0)   LB( I ) =0 

0021 

20 

CONTINUE 

0022 

PRINT  99 

0023 

PRINT  100,(IX(1,I),IX(2,K),IX(3,L),IX(4,M),IX(5,MM),I=1,3) 

0024 

PRI NT  1 01 , (LB( I ) , I =1 , 3) 

0025 

DO  40  1=1 , 127 

0026 

READ(9, 102, END=50)  N 

0027 

DO  25  J=1 , 3 

0028 

LB( J)=-2215, 925+XT( J)+IX(3, L) *5 , 39904+ I X( 4, M) *(-, 08788242) + 

*IX(5,MM)*1 17, 9473+N( 1 ) « 1 3 , 301 25+N ( 2 ) *340 , 7638+N ( 3 ) * ( -83 , 077 

*N(4)*368, 7023+N(5)*34, 89601 +N ( 6) « ( -1 77 , 5357) +N ( 7 ) *334 , 5793 

0029 

IF(LB( J) , LT, 0)  LB(J)=0 

0030 

25 

CONTINUE 

0031 

DO  30  J  =  1  ,  7 

0032 

IF(N( J) , EQ, n  A( J)=X 

0033 

30 

CONTINUE 

0034 

PRINT  103,  (A, LB( J) , J=1  ,  3) 

0035 

DO  35  J=1 , 7 

0036 

IF(A( J) ,EQ,X)  A(J)=B 

0037 

35 

CONTINUE 

0038 

40 

CONTINUE 

0039 

REWIND  9 

0040 

45 

CONTI NUE 

0041 

50 

STOP 

0042 

99 

FORMAT( ' 1 ' , 52X, 'M  I   N  E  D      AREA  S'/'  ',41X, 

*'NATIVE      SPECIES  PRODUCTION') 

0043 

100 

FORMAT( '0' , 3(35( ' -' ) , 9X)/'    ',3(':   PR      GS      K' , 

*'         NA      PH                   :',9X)/'    ',3(';    ' , 1 2, 2X, 4 ( 1 3, 2X) , 

*2X, ' ;           ;',9X)/'    ' , 3( ' ; ' , 27X, ' :   PR0D;',9X)/'  ', 

*3( ':', 7( '---+'), '     PER:',9X)/'    ',3(':TIL     SMTPS  ', 

*'FER  IRR  MUL    ST;   ACRE;',9X)/'    ', 3 ( 8 ('+---'),'--+' , 

*9X)  ) 

0044 

101 

FORMATC    '  ,  3(  '  ;  '  ,  27X,  16,  '  ;  '  ,  9X)  ) 

0045 

102 

F0RMAT(7I 1 ) 

0046 

103 

FORMAT ( '    '  , 3( ' ; ' , IX, Al , 6(3X, AI ) , IX, 16, ' ; '  ,  9X) ) 

0047 

END 
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(b)  Introduced  vegetation 


0001 
0002 
0003 

0004 
0005 
0006 
0007 
0008 
0009 
0010 
001  1 
0012 
0013 
0014 
0015 
0016 
0017 
0018 

0019 

0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 
0028 


0029 
0030 
0031 
0032 
0033 
0034 
0035 
0036 
0037 
0038 
0039 
0040 
0041 

0042 

0043 


0044 
0045 
0046 
0047 


10 

15 

»  X  « 


20 


DIMENSION  IX(5, 3) , LB(3) , A(7)  ,XT(3) 
INTEGER«2  N(7) 

DATA  IX/5, 50, 0, 0, 4, 1 5, 85 , 200 , 1 50, S , 25, 120,400,300,8/, 
«N/7«0/,X/'X'/, A/7*'    '/,B/'  '/ 
CALL  ERRSET  (208,256,-1,0) 
DO  45  MM=1 , 3 
DO  45  M=l , 3 
DO  45  L=l,3 
DO  45  K=1 , 3 
DO  20  1=1,3 
Z=IX(2,K) 

YPPR=(EXP(-(ABS( (5.0/7,0-1 .0)/.9)«*4,6))) 
IF( IX(2,K) .GT. 85)  GO  TO  10 

YPG1  =2510*(EXP(-(ABS(  (Z/85,  -1  .  )/.  1  6)  »«<4 )  )  ) +940 
XT(  I  )=YPPR«<YPG1  *  .  0061  8959«FL0AT (  I X ( 1 ,  I  )  )**1  .  6 
GO  TO  15 

YPG2=1200*(EXP(-(ABS( (Z/85. -1 . )/. 1 2) ««4) ) ) +2250 
XT(  I  )=YPPR«YPG2«  .  0061  8959'^FL0AT(  I  X(  1  ,  I)  )««1  .  6 
XT( I )=XT( . 17448 
FOR  ZERO  COMPUTATION.  «** 

LB( I )=-1 1 80. 820+XT( I )+IX(3, L) *4 . 08251 +IX( 4, M) * ( -1 .361515)+ 
«IX(5,MM)«143. 3554 
IF(LB( I ) . LT. 0)  LB(I)=0 
CONTINUE 
PRINT  99 

PRINT  100,  (IXd,  I)  ,  IX(2,K),  IX(3,L),  IX(4,M),  IX(5,MM)  ,  1  =  1,3) 
PRINT  101 , (LB(! ) , I =1 , 3) 
DO  40  I =1 , 127 
READO,  102,END  =  50)  N 
DO  25  J=1 , 3 

LB(J)=-1 180. 820+XT( J)+IX(3, L)*4. 08251 +IX(4,M)*(-1 .361515)+ 
«IX(5,MM)«143.  3554  +  N(  1  )  «  ( -307 .  9574 ) +N(  2)  «  ( -782 .  8676 )  ■«  N(  3)  >«99  .  46091  + 


25 


30 


35 
40 

45 
50 

99 

100 


5873+N(5)*418, 
LT.O)  LB(J)=0 


*N(4)*429 
IF(LB( J) 
CONTINUE 
DO  30  J=1,7 
IF(N( J) . EQ. 1 )  A( J)=X 
CONTINUE 

PRINT  103, (A, LB( J) , J=1 , 3) 

DO  35  J=1 , 7 

IF(A(J) .EQ.X)  A(J)=B 

CONTINUE 

CONTINUE 

REWIND  9 

CONTINUE 

STOP 

DEBUG  INIT(LB,XT) ,UNIT(2) 
FORMAT ( ' 1 ' , 52X, 'M  I  NED 
^'INTRODUCED  SP 


1331 +N( 6) « 159. 951 1 +N ( 7 ) « ( -78 . 4929 ) 


FORMAT ( ' 0 ' , 3 ( 35 ( ' - ' ) ,  9X ) / ' 
«'         NA      PH  : 
*2X, ' :  :',9X)/'  ',3( 


AREA  S'/'  ' 
E  C  I   E  S  PR 
PR  GS 
3(  '  : 
PROD 


C  T  I  0  N'  ) 


, 7( ' ---+' ) , 
IRR  MUL  ST 


*3(  '  : 
*'FER 
^  9X ) ) 

101  FORMATC  ',3( 

102  F0RMAT(7I1) 

103  FORMATC  ',3( 
END 


PER:  ' 
ACRE: ' 

27X, 16, 


.3( 
9X)/' 
'.27X, 
9X)/' 
9X)/' 


33X, 
0  D  U 
K', 

12,  2X,  4(  13,  2X) , 
, 9X ) / '    '  , 


3( ' : TIL     SM  TPS 
3(8( '+---' )^ '--+ 


9X)  ) 


IX, A1 ,6(3X, Al ) , IX, 16, ' : ' ,9X) ) 
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b.  Plant  cover  density  model 

(1)  Interactive  component  equation 


(Equation  4) 


Cover  =;  YP  * 


PR 
26 


-  1 


1  -  I 


N 


N 


-  e 


1  -  e 


Kind  of 
vegetation 

1.07686  (Native) 
0.98256  (introduced) 


Where 


N  =  1  +  NYP  *  (1.1397  *  e 


AGE 

10 


-  1 


0.88 


-  0.1397) 


(180  -  GS) 
108 


-  1 


0.52 
NYP  =  1.8  +  e 

I  =  0.38  +  (IYP/10)  *  AGE 


lYP  =  0.29  *  e 


GS 
180 


-  1 


0.5105 


-  0.2 


YP  =  YPGS  * 


(AGE  +  1) 
11 


-  1 


(1  -  IGS) 


_  e   {(1  -\gS)} 


1  -  e 


{(1  -\gs)} 


YPGS  =  100  *  e 


GS 
180 


-  1 


0.78 


IGS  =  0.1  +  0.23  *  e 


(180  -  GS) 
180 


0.36 


LIMITS 


5  <  PR    <  25  PR 
50  <  GS    <  180  GS 
0  <  AGE  <  7  AGE 


annual  precipitation,  inches 
growing  season,  days 
age  of  planting,  years 
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(2)  Additive  component  table 


(Table  2) 


Additive  plant  cover  density 
Component   (percent)  


Native  Introduced 


Soil  potassium  (K) 

^0.0892 

^0.0205 

Soil  sodium  (Na) 

1  

0.0186 

^0.0093 

Soil  pH 

2  

4.2 

7.6 

Tilliage  (Til) 

3.9 

10.0 

Seeding  method  (SM) 

broadcast 

0 

0 

drilled 

_ 

1.8 

+ 

1.0 

Topsoil  (TPS) 

+ 

6.0 

+ 

0.4 

Fertilizer  (FER) 

+ 

11.2 

+ 

30.9 

Irrigation  (IRR) 

+ 

0.6 

+ 

15.0 

Mulching  (MUL) 

+ 

15.2 

+ 

4.8 

Seeding  time  (ST) 

spring 

0 

0 

fall 

5.6 

+ 

19.0 

Intercept 

12.9 

3  _ 

74.3 

'Multiply  these  density  (percent)  values  by  parts  per  million  of  soil  com- 
ponent (potassium  or  sodium). 

^Multiply  these  density  values  by  pH  units  of  soil  acidity. 

^These  intercept  values  must  be  added  to  the  total  value  of  cover  density. 


LIMITS 

0  <  K    <  450 
0  <  Na  <  1000 
4  <  pH  <  9 
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(3)  FORTRAN  IV  computer  programs  for  estimating  plant  cover  density  (equation  4)  for: 
(a)  Native  vegetation 


0001 
0002 
0003 

0004 
0005 
0006 
0007 
0008 
0009 
001  0 
001  1 
0012 
0013 
001-1 
0015 
0016 
0017 
001  8 
0019 
0020 
0021 
0022 
0023 
0024 

0025 

0026 
0027 
0028 
0029 
0030 
0031 
0032 
0033 
0034 
0035 


0036 
0037 
0038 
0039 
0040 
0041 
0042 
0043 
0044 
0045 
0046 
0047 
0048 
0049 

0050 


DIMENSION  IX(5, 3) , LB(3) , A(7)  ,XT(3) 
INTEGER>k2  N(7) 

DATA  I X/5, 50, 0,0, 4, 15,85,200, 150,6,25, 120,400,300,8/, 
*N/7*0/,X/'X'/, A/7*'    '/,B/'  '/ 
CALL  ERRSET  (208,256,-1,0) 


DO 
DO 
DO 
DO 
DO 
V  = 


45 
45 
45 
45 
20 
IX(  1 


MM  = 
N=1 
L=1 
K=1 
1=1 
.,  I  ) 


Z=IX(2,K) 

AYP=1 . 8+EXP(-(ABS( ( ( 1 80-Z)/1 80, -1 


)/. 52)«*15) ) 


1 

1397)+1 


C  *  «  K 
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BYP=  .  29«  (EXP(  -(ABS(  (  Z/ 1  80  .  - 1  .  ) /  .  51  05 )  « «  1  2 )  )  ) -■ ,  2 
YPGS=100. »(EXP(-(ABS( ( Z/1 80 , - 1 , ) / . 78 ) * x 8 ) ) ) 
BGS=  ,  23*(EXP(  -(ABS(  (  (  1  80  .  ~Z ) /I  80  .  - 1  .  ) /  .  36 )  « >^6  ,  5 )  )  )  + 
AX  =  AYP*(  1  .  1397*(EXP(-(ABS(  (  5  ,  /  I  0  .  -1  .  )/ .  88  )  »>-5  .  8  )  )  )  ■- 
BX=(BYP/10. )*5, +.38 

CN=EXP( -(ABS( ( (5. +1 )/1 1 , -1 , )/( 1 . -BGS) )«*10) ) 

RN=EXP( -( ( 1 . /( 1 . -BGS) ) «*1 0) ) 

YP= ( (CN-RN)/( 1 . -RN) ) *YPGS 

TN  =  EXP(--(ABS(  (V/26,  -1  .  )/(  1  .  ~BX)  )**AX)  ) 

UN=EXP(-( ( 1 ./( 1 . -BX) )**AX) ) 

COV=( (TN-UN)/( 1 , -UN) )*YP 

XT( I )=CaV«1 . 07686 

FOR  ZERO  COMPUTATION.  ««* 

LB( I )  =  1 2. 8S572+XT( I)  +  IX(3,L)*. 0891  8672+ I  X ( 4 , M ) * ( - . 01 861335)  + 
*IX(5,MM)*(-4. 176856) 
IF(LB( I ) .LT.O)   LB( I )=0 
IF(LB( I ) .GT. 100)  LB( I )=100 
CONTINUE 


PRINT 
PRI  NT 
PRINT 
DO  40 
READ (9 


99 
100, 
101  , 

1=1  . 

,  1  02 


( IX( 1 , I ) 
(LB( I ) , I 
127 

, END=50) 


IX(2 
^  ,  3) 

N 


K),IX(3,L),IX(4,M),IX(5, MM), 1=1,3) 


DO  25  J=1 , 3 

LB( J)  =  1 2. 86572+XT( J)  +  IX(3, L) * , 089 1 8572+ I  X ( 4 , M ) * (  -  . 01 861335)  + 
*IX(5,MM)«(-4. 17685G)+N( 1 ) * ( -3 . 9291 49) +N( 2) * ( -1 . 8361 47) +N( 3) * 
*6.  056n78  +  N(4)*1  1  ,  1  9705-!-N  (  5 )  *  .  5945488  +  N  (  6 )    1  5  ,  !6765  +  N(7)* 
*(-5. 601713) 

0)  LB(J)=0 
100)  LB(J)=100 


LT 
GT 


25 


30 


35 
40 

45 
50 

99 


IF(LB( J) 
IF(LB( J) 
CONTI NUE 
DO  30  J=1 , 7 
IF(N( J) , EQ 
CONTI NUE 
PRI NT  1 03, (A, 
DO  35  J=1  ,  7 
IF(A( J) . EQ.X) 
CONTINUE 
CONTINUE 
REWIND  9 
CONTINUE 
STOP 

DEBUG  INIT(LB,XT) 
FORMAT ( ' 1 ' , 52X, 'M 
('NATIVE  S 


1 )  A( J)=X 


LB( J) , J=1 , 3) 


A( J)=B 


, UNIT (2) 
I   N  E  D 
P  E  C  I  E 


R  E  A  S'/'    ' , 41X, 
COVER') 


0051  100     FORMAT('0',3(35('-'),9X)/'    ',3(';   PR      GS      K' , 

*'         NA      PH         ;  ;',9X)/'    ',3(':    '  ,  I  2, 2X, 4 ( I  3, 2X) , 

*2X,':  ;',9X)/'    ' , 3( ' ; ' , 27X, ' ;  ;',9X)/'  ', 

*3( ';', 7( '---+'), '  %  ;',9X)/'  ',3(':TIL  SM  TPS  ', 
*'FER  IRR  MUL  ST ; COVER :', 9X )/ '  ', 3( 8 ('+--"'),'--+' , 
*9X)) 

0052  101     FORMATC    '  ,  3 (  '  :  '  ,  27X  ,  I  6  ,  '  ;  '  ,  9X )  ) 

0053  102  F0RMAT(7I1) 

0054  103     FORMATC    '  ,  3 (  '  ;  '  ,  1 X,  A1  ,  6 (  3X,  A1  ) ,  1  X  ,  1  6  ,  '  ;  '  ,  9X )  ) 

0055  END 
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(b)  Introduced  vegetation 


0001 
0002 
0003 

0004 
0005 
0006 
0007 
0008 
0009 
0010 
001  1 
0012 
0013 
0014 
0015 
0016 
0017 
001  8 
0019 
0020 
0021 
0022 
0023 
0024 

0025 

0026 
0027 
0028 
0029 
0030 
0031 
0032 
0033 
0034 
0035 


0036 
0037 
0038 
0039 
0040 
0041 
0042 
0043 
0044 
0045 
0046 
0047 
0048 
0049 
0050 

0051 


0052 
0053 
0054 
0055 


DIMENSION  IX(5, 3) , LB(3) , A(7)  ,XT(3) 
INTEGER«2  N(7) 

DATA  IX/5, 50, 0, 0, 4, 15, 85, 200, 150, 6, 25, 120, 400, 300, 8/ 
xN/7*0/,X/'X'/, A/7*'    '/,B/'  '/ 
CALL  ERRSET  (208,256,-1,0) 
DO  45  m=1  ,  3 
M=1,3 
3 
3 
3 


DO 
DO 
DO 
DO 


45 
45 
45 
20 


L=1 
K=1 
1  =  1 
V=IX( 1 , I) 
Z=IX(2,K) 
AYP=1 . 8+EXP( 


(ABS( ( ( 1 80-Z)/1 80, -1 .)/. 52)^*15) ) 


C  *  «  * 


20 


BYP=.29«(EXP(-(ABS( (Z/180. -1 .)/.5105)««12)))-, 
YPGS=100. *(EXP(-(ABS( (Z/1 80. -1 . ) / , 78 ) * « 8 ) ) ) 
BGS=.23«(EXP(-(ABS( ( ( 1 80 . -Z ) /I  80 . - 1  . ) / . 36 ) « *6 , 
AX  =  AYP«(  1  .  1397=^(EXP(-(ABS(  (5,/10,  -1  .  )/.  88)«*5, 
BX=(BYP/10.  )>^5.  +  .  38 

CN=EXP(-(ABS( ( (5. +1 )/1 1 . -1 . )/( 1 . -BGS) )*«10) ) 
RN=EXP(-( ( 1 ./( 1 . -BGS) )«*10) ) 
YP=( (CN-RN)/( 1 . -RN) )«YPGS 
TN=EXP(-(ABS( (V/26. -1 . )/( 1 
UN=EXP(-((1 ./(I .-BX))*«AX)) 
COV= ( ( TN-UN ) / ( 1 . -UN ) ) * YP 
XT( I )=COV« . 98256 
FOR  ZERO  COMPUTATION.  *** 
LB( I )=-74. 25282+XT( I )+IX(3 
*IX(5,MM)«7.S43977 
IF(LB( I ) . LT. 0)   LB( I )=0 
IF(LB( I) .GT. 100)   LB( I )=100 
CONTINUE 
PRINT  99 
PRINT 
PRINT 
DO  40 
READ (9 


5)  )  )  + 
8)  )  )- 


1 

1397)+1 


-BX) )**AX) ) 


L ) * . 02049003+ I X ( 4 , M ) * . 009275586  + 


100,  ( IX( 1 

101 .  (LB ( I 
I =1 , 127 
, 102, END=50) 

3 


, I ) , IX(2 
)  ,  I  =  1 , 3) 


K) , I X ( 3 , L ) , I X ( 4 , M ) , I X ( 5 , MM ) ,  I  =  1 , 3 ) 


N 


DO  25  J=1 

LB(  J)=-74.  25282+XT(  J)  +  IX(3,  D*  .  02049003+ 1  X(  4,  M)  >^ .  009275586  + 
*IX(5,MM)*7.  643977  +  N(  1  ) (  - 1  0  .  02594  ) +N  (  2 )  *  1  .  02138  +  N(3)* 

3622047+N(4)*30. 88224+N(5)>Kl5. 06255+N(6)*4. 778853+N(7)«18.96781 

IF(LB( J) . LT. 0)  LB(J)=0 

IF(LB( J) .GT. 100)  LB( J)=100 
25  CONTINUE 

DO  30  J=1 , 7 

IF(N( J) . EQ. 1 )  A( J)=X 
30  CONTINUE 

PRI NT  1 03. (A, LB( J) , J=1 , 3) 

DO  35  J=1 , 7 


IF(A( J) . EQ 
35  CONTINUE 
40  CONTINUE 
REWIND  9 
45  CONTINUE 
50  STOP 
99  FORMAT ( ' 1 ' 
«  '  I   N  T  R  0 
100     FORMAT  CO' 
*'  NA 


X)  A(J)=B 


, 52X, 'MINED      AREA  S 
DUCED  SPECIES 
, 3(35( ' -' ) , 9X)/'    '.3(':  PR 
PH         ;  : ' , 9X)/'    ' , 3( ' : 


/'  ',33X, 
COVER') 
GS      K'  , 

' , 12, 2X, 4( 13, 2X) 


*2X, ' :  :',9X)/'    ' , 3( ' ; ' , 27X, ' :  :',9X)/'  ', 

«3( ':', 7( '---+'), '     %     :',9X)/'    ',3(':TIL     SM  TPS 
«'FER  IRR  MUL     ST : COVER :', 9X )/ '       3( 8( '+---'),'--+' , 
*9X)  ) 

101  FORMATC    '  ,  3(  '  :  '  ,  27X,  16,  '  ;  '  ,  9X)  ) 

102  F0RMAT(7I1) 

103  FORMATC    '  ,  3C  C  ,  IX,  A1  ,  6(3X,  A1  )  ,  IX,  I  6,  '  C  ,  9X)  ) 
END 
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3.  Statistics  pertinent  to  each  model 


(Table  3) 


Statistics 


Prediction 
model 


Number  Coefficient 
of  of 

transects  determination 
(n)  ^R^) 


Standard 
error  of 
estimate 

(^y-x) 


Forage  production  83 
weight  on  unmined 
areas 

Plant  cover  density  83 
on  unmined  areas 


0.545 


0.785 


356  lbs/acre 


13.1  percent 


Forage  production 
weight  on  mined 
and  revegetated 
areas 

Native  vegeta- 
tion 44 
Introduced 
vegetation  33 

Plant  cover  den- 
sity on  mined  and 
revegetated  areas 
Native  vegeta- 
tion 44 
Introduced 
vegetation  33 


0.685 
0.735 


0.666 
0.615 


380  lbs/acre 
547  lbs/acre 


17.5  percent 
19.2  percent 


Packer,  Paul  E.,  Chester  E.  Jensen,  Edward  L.  Noble,  and  John  A.  Marshall.  1982. 
Models  to  estimate  revegetation  potentials  of  land  surface  mined  for  coal  in  the 
West.  USDA  For.  Serv.  Gen.  Tech.  Rep.  iNT-123,  25  p.  Intermt.  For.  and  Range 
Exp.  Stn.,  Ogden,  Utah  84401. 

Practical  criteria  for  measuring  success  of  revegetation  of  land  surface  mined  for 
coal  in  the  West,  namely,  the  amount  of  forage  produced  and  the  density  of  plant 
cover  developed,  are  affected  significantly  by  at  least  two  major  climatic  factors 
that  are  not  readily  susceptible  to  alteration  (precipitation  and  growing  season 
length);  by  three  properties  of  spoil  materials  that  are  subject  to  limited  modifica- 
tion through  management  (potassium,  sodium,  and  pH).  Revegetation  success  is 
also  influenced  by  seven  revegetation  treatments,  each  of  which  provides  at  least 
two  management  alternatives,  and  by  the  age  of  the  vegetation.  These  charac- 
teristics account  for  from  one-half  to  three-fourths  of  the  total  variance  in  forage 
production  and  plant  cover  density  on  these  revegetated  lands. 


KEYWORDS:  strip  mining,  surface  mining,  reclamation,  revegetation,  production 
 potentials,  cover  potentials  
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